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ABSTRACT 

We report on the discovery of the shortest period binary comprising a hot subdwarf star (CD— 30 
11223, GALEX J1411— 3053) and a massive unseen companion. Photometric data from the All Sky 
Automated Survey show ellipsoidal variations of the hot subdwarf primary and spectroscopic series 
revealed an orbital period of 70.5 minutes. The large velocity amplitude suggests the presence of a 
massive white dwarf in the system (M2/M0 > 0.77) assuming a canonical mass for the hot subdwarf 
(0.48 Mq), although a white dwarf mass as low as 0.75 Mq is allowable by postulating a subdwarf 
mass as low as 0.44 Mq. The amplitude of ellipsoidal variations and a high rotation velocity imposed 
a high- inclination to the system (i > 68°) and, possibly, observable secondary transits (i > 74°). At 
the lowest permissible inclination and assuming a subdwarf mass of ~ 0.48 Mq, the total mass of 
the system reaches the Chandrasekhar mass limit at 1.35 Mq and would exceed it for a subdwarf 
mass above 0.48 Mq. The system should be considered, like its sibling KPD 1930+2752, a candidate 
progenitor for a type la supernova. The system should become semi-detached and initiate mass 
transfer within « 30 Myr. 

Subject headings: binaries: close — white dwarfs 



1. INTRODUCTION 

A possible scenario for triggering type la supernovae 
involves a degenerate core collapse following th e merger 
of two white dwarfs (see a recent summa ry in iFerraricj 

mand recent population syntheses by iToonen et al.1 
) . A search for double degenerate stars that could 
progress toward such an outcome remains unfulfilled with 
mo st binaries either too lea n or too wide (SPY program, 
see iNapiwotzki et al.l 120011), with the possibl e exception 
of WD 2020-425 (|Napiwotzki et al. 1 120051) , A more 
promising variation on this scenario involves the inter- 
action of a hot subdwarf star with a close white dwarf 
companion. Mass transfered from a He-star onto a CO 
white dwarf may ignite and even lead to detonation of 
the carbon core. At an earl y stage these ob j ects ap- 
pear as AM CVn systems (see lYunge lson 2008; S olheiml 
2010). Hot subdwarfs, also referred to as extreme hori- 
zontal branch (EHB) stars are core helium burning stars 
that formed prefere ntially in the aftermath of close, post- 
interacting binaries (IHan et al.ll2002t [20031. The identifi- 
cation of the short-peri od subdwarf plus wh i te dwarf bi- 
nary KPD 1930+2752 (jMaxted et al.l I200CH : IGeier et al.l 
I2007h lent support to the idea, and the expectation of 
conservative mass transfer gave KPD 1930+2752 the 
character of a proper type la candidate progenitor. A 
search for other massive subluminous companions to hot 
subdwarf stars has so far failed to pr ovide further can- 
didate s (MUCHFUSS program, see IGeier et al.l f201ll 
1201231 . 

In this letter, we report on the discovery of a close 
companion to the hot subdwarf CD— 30 11223. The or- 



Based on observations made with ESO telescopes at the La 
Silla Paranal Observatory under programme ID 83.D-0540, 85. D- 
0866, 089.D-0864. 



bital period is the shortest known to date for such sys- 
tems. The observations were carried-out as part of an 
intermediate-dispersion survey of ^40 ho t, ultraviolet- 
selecte d su bdwarf stars from th e catalog of iVennes et al.1 
(|2011f ) and lNemeth et all (|2012l ) and aimed at constrain- 
ing the fract ion of close, e volve d binaries in the EHB 
population. iKawka et al.1 (|2010l ) described the orbital 
properties of the first two systems discovered as part of 
this survey. 

The hot hydrogen-rich subdwarf (sdB) CD-30 11223 is 
relatively bright (V « 12) and nearby (d = 250 - 400 pc) 
and was first classified by IVennes et al.l (|2011[) in spec- 
troscopic follow-up observations of ultraviolet sources 
detected in the Galaxy Evolution Explorer (GALEX). 
Based on low-dispersion spectra, they originally deter- 
mined T eS = 26500 ± 2700 K, \ogg = 5.36 ± 0.26, and an 
abundance of helium log He /H = —1.34 + _ 0.24. These 
parameters were refined by iNemeth et al.l (|2012f ) using 
higher-dispersion spectra: T e ff= 3015Ql 310 K, logg = 

5.72lg$, logHe/H = -1.58+^. INemeth et all (l2012Tl 
noted that temperature and surface gravity measure- 
ments in low dispersion spectra were systematically off- 
set from measurements obtained at hi gher-dispersion 
and a re probably less reliable. Moreover, INemeth et al.l 
(2012) reported the detection of nitrogen lines and mea- 
sured an abundance logN/H = —4.40^0 67- They found 
no evidence of carbon or oxygen in the atmosphere 
of this star. Remarkably, large radial velocity excur- 
sions between the two series of spectra were evident and 
prompted further investigations. Independently, photo- 
metric variations in a large ASAS data set confirmed our 
suspicions that this object is in a close, evolved binary. 
The nature of the variations (reflection effect or tidal dis- 
tortion) as well as the nature of the companion may be 
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established with joint phasing of photometric and spec- 
troscopic data. 

In the following sections we describe the spectroscopic 
and photometric data (Section 2) and their analysis (Sec- 
tion 3). Section 3.1 gives the results of a frequency anal- 
ysis of the photometric data and Section 3.2 provides 
corroborating evidence of the short orbital period based 
on spectroscopic data. We also measured the mass func- 
tion of the secondary and rotation velocity of the primary 
and analyzed the ellipsoidal variations to constrain the 
nature of the secondary star (Section 3.3). Finally, we 
summarize in Section 4. 

2. OBSERVATIONS 

The first identification spectra of the hot subdwarf 
CD— 30 11223 were obtained with the European South- 
ern Observatory (ESO) Faint Object Spectrograph and 
Camera (EFOSC2) attached to the 3.6 m New Technol- 
ogy Telescope (NTT) at La Silla Observatory on UT 2009 
August 24 (ESO P83 program). We used grism no. 11 
which has 300 lines per mm and a blaze wavelength at 
4000 A. The slit width was set to 1 arcsec, which re- 
sulted in a spectral resolution of ~14 A and a coverage 
from 3700 to 7250 A. We re-observed CD-30 11223 at 
the NTT on UT 2010 September 21 (ESO P85 program). 
We employed grism no. 7 (600 lines per mm) with a dis- 
persion of 1.96 A per pixel. The slit width was set at 
1 arcsec resulting in a resolution of ^8 A . Although 
the 2009 and 2010 spectra only offered modest velocity 
accuracy (50 to 60 kms -1 ) the difference in measure- 
ments was close to 800 kms -1 suggesting the presence of 
a massive, close companion to the hot subdwarf. Con- 
sequently, we secured a set of intermediate dispersion 
spectra, first at Siding Spring Observatory (SSO) using 
the 2.3m telescope (April 2012) and at ESO using the 
NTT (August-September 2012). 

We obtained a set of three spectra on UT 2012 April 
29 and 30 with th e Wide Field Spectrograph (WiFeS 
IDopita et al.l H007) attached to the 2.3m telescope at 
SSO. We used the B3000 and R7000 gratings which pro- 
vided spectral ranges of 3200-5900 A at a resolution 
R = A/AA = 3000, and 5300-7000 A at R = 7000, re- 
spectively. We used the RT560 dichroic beam-splitter to 
separate the incoming light into its red and blue compo- 
nents. We maximized the signal-to-noise of each observa- 
tion by extracting the spectrum from the most significant 
(< 6) traces. Each trace was wavelength and flux cali- 
brated prior to co-addition. The spectra were wavelength 
calibrated using NeAr arc spectra that were obtained fol- 
lowing each observation. The exposure time was 600 s. 
The two exposures acquired on April 29 were separated 
by only 45 minutes but revealed a velocity shift close to 
500 kms" 1 . 

Finally, we obtained complete orbital cycles at the 
NTT from UT 2012 August 30 to September 2 (ESO 
P89 program). We used grism no. 20 centered on Ha 
and with a dispersion of 0.55 A per pixel in the spec- 
tral range from 6040 to 7140 A. We set the slit width 
to 0.7 arcsecond resulting in a ~ 2 A resolution. We set 
the exposure times to 120 s in order to minimize orbital 
smearing. We obtained 8 spectra on the night of August 
31 and 6 spectra on September 1. The spectra revealed 
short (« 70 min) high- amplitude velocity variations and 
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Fig. 1. — (Bottom panel) Frequency analysis of the ASAS photo- 
metric data (inverse of the normalized reduced \r versus frequency) 
revealing one significant frequency near 40.83 cycles d" 1 (full line). 
The 1-cr significance level (horizontal dashed line) sets a formal sta- 
tistical error on the chosen frequency. A weaker signal (not shown) 
is also present at / p hot/2. (Upper panel) ASAS photometric data 
folded on the best frequency and best fit at / p hot/2- 

appeared to cover several orbits. 

The field near CD— 30 11223 was a lso covered with 
the A ll Sky Automated Survey (ASAS, IPoimanskilli~997L 
l2001f ) on numerous occasions. A total of 1060 images 
were obtained as part of ASAS using the photometric 
V band. The star CD— 30 11223 was not immediately 
identified as a variable star. We retrieved the photo- 
metric measurements using the ASAS all star catalogue 
pageQ. 

3. ANALYSIS 
3.1. Photometric frequency analysis 

The ASAS measurements span nearly 3200 days from 
UT 2000 December 24 to 2009 September 13. Figure [Q 
shows a frequency analysis of the ASAS data. Fitting 
a sinusoidal function to the data set unambiguously se- 
lected a signal frequency / phot = 40.83377 ± 0.00003d" 1 
with a semi-amplitude Am/2 = 0.047 ± 0.004 mag and 
mean magnitude my — 12.342±0.003. A reflection effect 
on a fainter companion would display sinusoidal varia- 
tions on the orbital period P orD (= 2/r^ t ), while tidal 
distortion of the primary star would induce sinusoidal 
variations on half the orbital period P rb/2 (= /phot) as 
shown in Figure [TJ An analysis of radial velocity data 
should help elucidate the nature of these variations. 

3.2. Radial velocities and secondary mass function 

All spectra were obtained at high signal-to-noise ra- 
tio and revealed well defined B aimer and Hei lines. We 
fitted the Ha (or H/3 in the P85 data) line cores using 
a Gaussian function. We estimated the errors at 50-60 

2 http:/ /www.astrouw.edu.pl/asas/?page=main 
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kms 1 in the low-dispersion data (P83 and P85) and 5- 
10 kms -1 in the intermediate dispersion data (P89 and 
SSO). Table[T]lists all radial velocity measurements. The 
heliocentric velocity is given as a function of the HJD 
time at mid-exposure. Folding the data on half the pho- 
tometric frequency proved successful showing that the 
variations are due to tidal distortion of the bright hot 
subdwarf primary. The corresponding ephemeris is 

P = 0.04897906 ± 0.00000004 d 

T = 1967.8426 ±0.0012 (1) 

where To is the time of passage of the hot subdwarf pri- 
mary at inferior conjunction. 

Figure [5] shows the radial velocity data folded on the 
(above) ephemeris. The ASAS data preceded most of 
our radial velocity data by ~ 1000 days and the accu- 
mulated error in the predicted phase is ~ 1.7 minutes. 
In fact, our radial velocity measurements are offset from 
the predicted radial velocity curve by only 0.0246 phase 
(or ~ 1.7 minutes) and firmly link the photometric vari- 
ations to tidal distortion of the subdwarf primary. The 
residual in the more precise P89 velocity measurements 
is 12 kms -1 and is consistent with the estimated errors 
of 10 kms . The velocity semi- amplitude and mean 
velocity of the hot subdwarf are 
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= 31.5 ± 1.3 kms" 



#i(sdB) = 386.9 ± 1.9 kms , 

and the corresponding mass function for the unseen com- 
panion is 

/ 2 =0.294 ±0.004 M e . 

Finally, we generated a phase-corrected, co-added spec- 
trum of the Ha and HeiA6678 using all ESO P89 spec- 
tra. The exposure time limited orbital smearing to <50 
kms -1 , or nearly half of the spectral resolution of the 
P89 data (90 kms -1 ). A rotation velocity comparable 
or in excess of 100 kms" 1 should be readily measurable. 

3.3. Properties of the primary and nature of the 
secondary 

Approximating the full amp litude of the ellipsoidal 
variations by (see lWarnenll995[ ): 

Am = 3.5 5^ — sin 2 i 

where q = M1/M2, Mi and R\ are the mass and radius 
of the hot subdwarf primary, and M2 is the mass of the 
unseen companion. The subdwarf mass is constrained in- 
dependently but Af 2 and i are treated as free parameters. 
The semi-major axis a follows Kepler's third law. The 
gravity darkening and limb darke ning coe f ficient s were 
held fixed at values suggested by iWarnerl ()1995[ ). e.g., 
fj, g = 1 and /i = 0.36. 

The subdwarf radius may be obtained from Ri /R@ = 
y/Mx/ig/go). Adopting log q from lNemeth et al.l (|2012ft 
and assuming Mi = 0.48 M Q , we estimated the subdwarf 
radius Ri = 0.159l ' 07 Rq. We may also set an indepen- 
dent and more reliable constraint on the primary radius 
using a measurement of the rotation velocity (Fig [3]). 
Fitting synthetic Ha profiles to co-added P89 spectra 
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Fig. 2. — Radial velocity measurements (open squares: P83, 
P85 and SSO; full squares: P89) folded on the photometric-based 
ephemeris (Equation 1) compared to the predicted radial veloc- 
ity curve (full line). The velocity measurements are offset from 
the predicted phase by only 0.0246 (dashed line) consistent with 
accumulated errors in the photometric ephemeris. 
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Fig. 3. — Co-added, intermediate dispersion spectra (grey lines) 
and best-fitting, rotationally broadened synthetic spectra (full 
lines) near (bottom) Ha and (top) HelA6678. The synthetic spec- 
tra without rotation broadening are shown with dashed lines. The 
effect of a projected rotation velocity of 164 kms -1 is evident in 
the line cores. 



and adopting best-fit model parameters (|Nemeth et al.l 
2012) we measured v ro t sinz = 164 ±5 kms 1 . The mod- 
els were first convolved with a Gaussian function with a 
full-width half-maximum set at 2A. Also, the effect of 
orbital smearing was taken into account by smoothing 
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TABLE 1 
Radial velocities. 
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ESO NTT (P83) 


5067.49158 
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ESO NTT (P89) 


6171.46397 


-296.2 


ESO NTT (P89) 


6171.54961 


-142.9 


ESO NTT (P83) 


5067.49268 


436.0 


ESO NTT (P89) 


6171.47649 


225.6 


ESO NTT (P89) 


6172.47112 


283.2 


ESO NTT (P85) 


5460.48093 


-377.2 


ESO NTT (P89) 


6171.48926 


339.3 


ESO NTT (P89) 


6172.48510 


-329.0 


ESO NTT (P85) 


5460.48412 


-313.1 


ESO NTT (P89) 


6171.50260 


-245.3 


ESO NTT (P89) 


6172.49843 


-61.8 


SSO 2.3-m 


6047.25139 


-340.8 


ESO NTT (P89) 


6171.51655 


-196.1 


ESO NTT (P89) 


6172.51202 


415.7 


SSO 2.3-m 


6047.28282 


145.6 


ESO NTT (P89) 


6171.53024 


388.9 


ESO NTT (P89) 


6172.52588 


-25.8 


SSO 2.3-m 


6048.01768 


165.2 


ESO NTT (P89) 


6171.54290 


161.6 


ESO NTT (P89) 


6172.53945 


-334.2 



the models with a box-car f unctio n with a full-width of 
50 kms" 1 . IGeier fc Heberl (pOlU ) found that all subd- 
warfs from their sample, single or in binaries with a pe- 
riod longer than 1.2 d are slow rotators (u ro tsini < 10 
kms" 1 ). Clearly, CD-30 11223 is exceptional. We 
verified our solution by convolving with the rotation 
broadening function a template built using a subdwarf 
(GALEX J0206+1438) similar to CD-30 11223 and ob- 
served with the same set-up. We measured an identical 
value of v ro t sinz = 164 ± 15 kms -1 but with larger er- 
ror bars possibly because of a slight parameter mismatch 
between the two stars. Taking observed ellipsoidal vari- 
ations as evidence that the subdwarf is tidally locked to 
the primary, the rotation and orbital velocities are re- 
lated: 

Vrot sin i _ i?i 
Ki ai ' 

and since a = ai + a 2 — a\{\ + q) then 

a v ro t sin i 



where q = Afi/M 2 = a 2 /ai. The resulting pri- 
mary radius is consistent with the gravitational radius 
(0.159 Rq) for all probable inclinations (> 50°). 

Figure 0] shows a set of constraints placed on the free 
parameters M 2 and i. The mass function defines a nar- 
row constraint for a given subdwarf mass: The minimum 
secondary mass ranges from 0.75 to 0.77 Mq for a sub- 
dwarf mass ranging from 0.44 to 0.48 M©. Combining 
mass function and tidal constraints at M\ — 0.48 M©, 
the secondary mass ranges from 0.77 to 0.87 M© imply- 
ing a total mass of 1.25 — 1.35 M©. The effective tem- 
perature and surface gravity of CD— 30 11223 place it 
squarely amongst the bulk of sdB stars that are best de- 
scribed with the models of iDorman et al.l (|1993[ ) within 
a narrow mass range centered on ~ 0.48 M©. Moreover, 
tests showed that constraints placed by ellipsoidal vari- 
ations and by the mass function are incompatible below 
Mi/Mq » 0.44. Therefore, assuming a mass close to 
0.48 M© is probably justified and the possibility that past 
interactions may have left an extremely low mass rem- 
nant (« 0.3 Mq, see, e.g.. IHeber et all 120031 ) does not 
need to be considered. 

Uncertainties on the rotation velocity measurements 
are possibly affected by systematic effects such as over- or 
under-estimated spectral resolution or orbital smearing. 
Assuming a rotation velocity as low as 150 kms -1 would 
increase the upper mass limit to 1.05 M© and push the 
inclination as low as 57°. 

The possibility of shallow (Am » (R 2 /Ri) 2 < 5 
mmag) secondary transits cannot be ruled out and tran- 
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Fig. 4. — Constraints on the secondary mass (M2) as a function 
of the orbital inclination. First, M2 is obtained from the mea- 
sured mass function fa assuming Mi/Mq = 0.48 (thick full line). 
The range of primary mass is dictated by evolutionary consi dera- 
tions, such as the evolutionary t racks o f Dorman ct al. ( 1993) that 
generates remnants of 0.48 Mq (Zhang ct al. 2009; Fontaine ct al. 
12013) . Next, a relation between M2 and the inclination is derived 
from the measured amplitude of the ellipsoidal variations assuming 
Mi = 0.48 (thin lines from bottom to top Am =0.086, 0.094, 0.102 
magnitude and shaded area). The diagram is split between areas 
with detectable and non-detectable transits. 

sit timing would allow us to constrain accurately the in- 
clination and stellar radii. Transits may become observ- 
able at an inclination i such that tan (90° — i) < R\/a 
or i > 74° (where a s» 0.6 Rq). It is worth noting that 
all solutions provide a primary radius smaller than its 
Roche lobe radius (Rl ~ 0.2i?©) excluding the possibil- 
ity of present-day interaction. 

4. SUMMARY AND CONCLUSIONS 

We found that the hot subdwarf CD-30 11223 is in a 
close 70.5 minutes binary with an unseen companion pre- 
sumed to be a white dwarf. The orbital properties and 
noted photometric variations suggest the presence of a 
companion with a minimum mass of 0.75 M© assuming 
a minimum mass of 0.44 M© for the primary. However, 
the measured parameters of the sdB CD— 30 11223 favor 
a canonical mass close to 0.48 M© that would correspond 
to a minimum mass of 0.77 M© for the secondary. Con- 
straining the radius of the primary using a measurement 
of its rotation velocity and fitting the amplitude of the 
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ellipsoidal variations set an upper limit of 0.87 Mq for 
the secondary mass at an inclination i rs 68° and a total 
mass of 1.35 Mq for the system. However, systematic 
effects in the measurement of the rotation velocity and 
the possibility that the primary mass may exceed 0.48 
Mq imply that the system mass may exceed the Chan- 
drasekhar mass limit (Mch)- Observations of transits 
would help constrain the inclination of the system and 

the primary radius . 

Following IRitterl (|1986T) , the system will initiate mass 
transfer in 25 to 28 Myr with a contact period of 49 min- 
utes. While transferring mass, the system would take the 
appearance of an AM CVn sy stem that could eventually 
trigg er a type la supernova (jYungelsonl 120081 : iSolheiml 
120101) . With a total life-time on the EHB that must ex- 
ceed 100 Myr ()Dorman et al.lH993[ ) and the short time 
remaining before contact (< 30 Myr), it is likely that a 
substantial amount of helium may remain in the envelope 
of the EHB star that w ould then consti t ute a generous 
helium donor. Following Nelema ns et al.l ()2001[ ). the out- 
come depends on the mass of the accreting white dwarf 
and the amount of helium available from the donor. The 
minimum mass of the white dwarf (> 0.75 Mq) exceeds 
a threshold (0.6 Mq) below which the system would un- 
dergo nova-like, helium-flash episodes. Instead, should 
enough helium (> 0.3 Mq) accumulate on the massive 
white dwarf after the onset of stable mass transfer the 



system may suffer helium detonation that could light-up 
the C/O core triggering a type la supernova. On the 
other hand, it is possible that the EHB star could ex- 
haust most of its helium supply before contact and merge 
with its massive white dwarf companion. Should the to- 
tal mass of the system exceed Mch the remnant could 
also trigger a type la supernova. Alternatively, should 
the system mass prove insufficient, the immediate out- 
come would be a hot, ultramass i ve wh ite dwarf such as 
those sampled by IVennes et al.l (|1997f ) in the Extreme 
Ultraviolet Explo rer survey. 

We noted that IGeier et all (|2012bD independently re- 
ported on this system and reached similar conclusions. 

S.V. and A.K. acknowledge support from the Grant 
Agency of the Czech Republic (GA CR P209/10/0967, 
GA CR P209/12/0217). This work was also supported 
by the project RVO:67985815 in the Czech Republic. We 
thank the referee for helping us clarify the evolutionary 
prospects of the system. Also, we thank A. Williams 
(Perth Observatory) for his contribution to the project. 

Note added in proof: L. Yungelson informed us 
that, following Fink et al. (2010, A&A, 514, 53), a su- 
pernova may be triggered by an accreted mass of he- 
lium as low as 0.1 Mq for a 0.8 Mq white dwarf, thereby 
increasing the odds of such an event taking place in 
CD-30 11223. 
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